Directionality of THz emission from photoinduced gas plasmas 
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Forward and backward THz emission by ionizing two-color laser pulses in gas is investigated by means of a 
simple semi-analytical model based on Jefimenko's equation and rigorous Maxwell simulations in one and two 
dimensions. We find the emission in backward direction having a much smaller spectral bandwidth than in 
forward direction and explain this by interference effects. Forward THz radiation is generated predominantly 
at the ionization front and thus almost not affected by the opacity of the plasma, in excellent agreement with 
results obtained from a unidirectional pulse propagation model. © 2013 Optical Society of America 
OCIS codes: 260.5210, 300.6270, 320.7120, 350.5400 
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Generation of radiation in the THz range and control- 
ling its spectrum and direction of emission is crucial for 
applications reaching from nonlinear THz spectroscopy 
to biomedical and security imaging. Several methods for 
THz generation, such as photoconductive switches or op- 
tical rectification in second order nonlinear crystals are 
limited in achievable THz field amplitudes due to satu- 
ration or material damage for high input intensities. An 
alternative setup [1-6], where an ionizing two-color laser 
pulse is focused into a gas cell has attracted much in- 
terest, since the obtained THz pulses are characterized 
by comparably high amplitudes and broad spectra at 
the same time. The THz emission can be explained by a 
low frequency component in the plasma current caused 
by the asymmetric two-color laser field [2], where the 
stepwise modulation of the electron current by two-color 
tunneling ionization plays a crucial role [6]. Within this 
model, good agreement of forward (FW) emitted THz 
radiation in experiment and simulation was obtained re- 
cently [6]. However, to the best of our knowledge, an 
investigation of possible backward (BW) THz emission 
due to this mechanism is still missing. 

In this Letter we present a combined analytical and 
numerical study of THz generation by ionizing two-color 
femtosecond pulses in a gas. We reveal the governing ef- 
fect for differences in FW and BW emitted THz spectra 
to be interference of signals from different spatial posi- 
tions in the plasma channel. Furthermore, we find that 
opacity of the plasma is of minor influence on FW THz 
emission in focused pump geometries. 

Let us first develop a simple picture to explain the 
main difference between FW and BW emitted fields. A 
simplified though generic setup is shown in Fig. 1. We 
assume a plasma line-source of length L on the z-axis, 
which is meant to be produced by a laser pulse propa- 
gating in positive 2-direction (FW). At each point, the 
pump pulse ionizes the medium, the generated free elec- 
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Fig. 1. (a) Schematic illustration of the interference being 
responsible for shaping FW and BW emission. A plasma 
line-source of length L created by a propagating two- 
color pump pulse emits radiation. Blue circles centered 
around exemplary point-sources represent planes of con- 
stant phase of elementary spherical waves. In FW direc- 
tion (c), spherical waves interfere constructively for all 
wavelengths and the resulting on-axis spectrum is pro- 
portional to the single emitter spectrum. In BW direc- 
tion (b), the spectral form-factor Eq. (2) depletes wave- 
lengths smaller than the source length (here L = 30 /iiii). 



trons are accelerated in the laser field and thus build up a 
current J(t), which in turn emits electromagnetic radia- 
tion. Due to the propagation of the ionizing pump pulse, 
points with smaller z-coordinates emit radiation earlier 
than the ones with larger z-coordinates. If all propaga- 
tion effects (diffraction, dispersion, etc.) are neglected, 
the pulse moves unchanged and all constituent points 
along the plasma line emit the same field, but shifted in 
time. This is illustrated in Fig. 1(a), where a snapshot 
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of the emitting plasma line is shown. Each blue circle is 
centered around a different point-source, representing a 
plane of constant phase of an elementary spherical wave 
for one wavelength A. To record the resulting spectra on 
a screen located before (BW) or behind (FW) the line, 
Jcfimenkos's equation [7] 
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is used to calculate the electric field emitted from a given 
current distribution. Here, R = |r — f^j is the distance 
between the emitting point-source and the screen. 

In this simplified situation, the FW on-axis spectrum 
is just proportional to the spectrum of a single point 
emitter |£'g (A~^)|, because contributions from emitters 
along the plasma line add up constructively for all wave- 
lengths A [see Fig. 1(c)]. In contrast, the BW on-axis 
spectrum differs significantly from the single emitter 
spectrum [see Fig. 1(b)], since radiation from different 
emitters is superposed with different temporal delays. 
The influence of these delays on the BW spectrum can 
be accounted for by a wavelength dependent form-factor 

/(A-i) = Ei(^(z - L/2)) - Ei(^(z + i/2)), (2) 

which can be obtained directly from Eq. (1) by assum- 
ing a uniform current density being temporally shifted 

by z' /c along the line. Here, Ei(a;) ~ j-oo "^^/'^ dx' 
is the exponential integral and z is the position of the 
screen. The on-axis BW spectrum is then proportional 
to \f{\-^)E^^{\-^)\. Figure 1(b) reveals that the BW 
spectrum gets depleted for wavelengths A <^ L, whereas 
it coincides with the FW spectrum for \ ^ L, because 
for these wavelengths the plasma line appears as a point 
source. 

We now confront our predictions with rigorous two- 
dimensional (2D) Maxwell simulations. We consider two- 
color input pump pulses 
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with amplitude A, beam width w, pulse duration at, 
r± = X (2D) or r± = x"^ + y-^ (3D), relative strength 
of fundamental (wq ~ '^nvQ, vq = 375 THz) and second 
harmonic ^, relative phase ip, being focused into argon 
gas at atmospheric pressure. The density of generated 
free electrons obeys dtPe = W^st(^^) [pat — Pe{t)], where 
Pat is the neutral atomic density and Wst{E) a field 
dependent tunneling ionization rate [5,6]. Further on, 
we assume zero velocity for newly born electrons which 
are then accelerated in the electric field, leading to the 
current density [2,5,8] 
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We model 2D transversal electric (TE) field evolution 
by means of the finite difference time domain (FDTD) 
method [9], solving Maxwell's equations 
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Linear dispersive properties of argon are included via 
Dy{uj) = eon'^(uj)Ey(uj), where we use the refractive in- 
dex n{uj) given in [10]. In all FDTD calculations we 
neglect nonlinear contributions to the polarization, be- 
cause they play a minor role for THz generation in the 
present configuration. Figure 2 shows results for a fo- 
cused (/ = 125 pm) pump pulse with A = 13.6 GV/m, 
w = 32 pm, at = 34 fs, ^ = 0.2, and (p = 7r/2. In ex- 
cellent agreement with our theoretical predictions, the 
BW spectrum [Fig. 2(a)] is very narrow compared to 
the FW one [Fig. 2(b)]. Moreover, estimating the length 
L of the plasma channel to be 30 pm [Fig. 2(d)] and 
evaluating Eq. (2) leads to a BW spectral 1/e- width of 
Av = 5 THz, which coincides with the numerical ob- 
tained one at z = —60 pm. In time domain, the ratio 
between on-axis BW and FW THz field amplitude in 
Fig. 2(d) is about 0.1. To further validate the interference 
to be the governing effect shaping the BW spectrum, we 
linearly propagate a two-color Gaussian pump pulse in 
3D geometry and use Jefimenko's Eq. (1) to compute 
the BW spectrum (3D Jefimenko approach). In order to 
get a comparable plasma channel (length and free elec- 
tron density), wc adjust initial width and amplitude to 
w ^ 25 pm and A~ 6 GV/m compared to the previous 
2D configuration. The inset in Fig. 2(a) shows the re- 
sulting BW spectrum, reproducing again the predicted 
as well as simulated spectral width of 5 THz. 

In the simple picture we have developed in Fig. 1 and 
Eqs. (1) and (2) we completely ignore the opacity of the 
generated plasma for frequencies below the plasma fre- 
quency 2Tn'p = ■\/ Peq^/f-ome- In the 2D examples above 
we find w 45 THz, but the generated plasma channel 
is small in transverse direction [w 3 pm, see Fig. 2(d)] 
and THz radiation strongly diffracts, so wc find our spec- 
tra still governed by the predicted interference effects. 
However, this may change if larger plasma volumes are 
created. In order to check this issue we resort to the ex- 
treme case of one-dimensional (ID) configuration, where 
we neglect all transverse spatial dependencies. Two-color 
pump pulses are launched in vacuum and hit a 1 mm 
thick argon layer at z = mm. This setup enables us 
to record BW spectra in FDTD simulations. For addi- 
tional comparison, wc simulate the FW fields by using a 
unidirectional pulse propagation equation (UPPE) [11]. 
This approach describes pump as well as the emitted 
THz fields in FW direction and was shown to reproduce 
experimental results [6] . In our UPPE code, we can also 
include third order nonlinear polarization (Kerr effect), 
which is neglected in the FDTD simulations. It is worth 
noticing that the UPPE approach does not take into 
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Fig. 2. 2D FDTD simulations: THz spectra of (a) BW 
and (b) FW emitted fields. The inset shows the on-axis 
BW spectrum obtained from a linearly propagated Gaus- 
sian pump pulse and Eq. (1) (3D Jefimenko approach, see 
text), (c) On-axis BW spectral width vs. plasma channel 
length L from Eq. (2) (solid line), 2D FDTD calculations 
(blue circles), and 3D Jefimenko approach (red crosses), 
(d) Snapshot of the emitted THz fields (< 100 THz) 
and the plasma channel, illustrating the strong (weak) 
emission in FW (BW) direction. 
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Fig. 3. ID simulations: FW THz spectra for (a) A = 
31 GV/m at 1 bar and (b) A = 46 GV/m at 5 bar from 
FDTD simulations (red line) and UPPE simulations with 
(green line) and without (blue line) third order nonlin- 
ear polarization included. The insets show the complete 
spectra, (c) illustrates that THz emission (green line) 
takes place mainly at the ionization front (blue line) . (d) 
shows BW THz spectra corresponding to (a) (blue line) 
and (b) (green line) obtained from FDTD simulations. 
All spectra are normalized to \Ein [v = 375 THz)|. 



account the opacity of the plasma, whereas it is natu- 
rally included in the FDTD simulations. Figure 3 shows 
pump pulses with at = 34 fs, ^ = 0.2, and Lp = 7r/2 
and amplitudes A = 31 GV/m and 46 GV/m launched 
into argon at pressure p = 1 bar and p = 5 bar, re- 
spectively. The plasma, generated in the 1 mm long ar- 
gon layer, is opaque for frequencies below Vp = 22 THz 
{A = 31 GV/m, pif ^'^ = 6 X lO^^ m'^) and Vp = 104 THz 
[a = 46 GV/m, p^^'' = 13.5 x lO^^ m-^). it turns 
out that for the FW emitted THz radiation the opac- 
ity of the plasma plays almost no role, we obtain ex- 
cellent agreement between FDTD and UPPE simula- 
tions [see Fig. 3(a,b)]. The reason for this effect is de- 
tailed in Fig. 3(c): the THz field is emitted before the 
plasma builds up, thus no damping influences its further 
FW propagation. In contrast, in BW direction plasma- 
opacity is important and only frequencies above the cor- 
responding plasma frequencies can propagate through 
the already built up plasma. Therefore, in ID geometry 
Eq. (2) is applicable only for frequencies above v > Vp, 
and for a plasma length of 1 mm entirely destructive in- 
terference is expected. Indeed, this can be seen from the 
cutoffs of the BW spectra at the corresponding plasma 
frequencies in Fig. 3(d). Below these frequencies i^p, ra- 
diation originates from the front-side of the plasma layer 
only, since contributions from deeper inside are damped. 

In conclusion, we explained the differences in FW and 
BW emitted THz spectra from gas plasmas observed in 
full Maxwell numerical simulations by simple interfer- 



ence effects. Forward spectra are well described by uni- 
directional models and opacity of the plasma has negli- 
gible influence. We believe that the predicted sensitivity 
of the width of the BW THz spectrum on the length of 
the plasma channel could be tested experimentally in a 
straight-forward manner. 
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